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ABSTRACT 

NGC 6240 is a pair of colliding disk galaxies, each with a black hole in its core. We have used laser 
guide star adaptive optics on the Keck II telescope to obtain high-resolution 0.06") near-infrared 
integral-field spectra of the region surrounding the supermassive black hole in the south nucleus of this 
galaxy merger. We use the K-band CO absorption bandheads to trace stellar kinematics. We obtain 
a spatial resolution of about 20 pc and thus directly resolve the sphere of gravitational influence of the 
massive black hole. We explore two different methods to measure the black hole mass. Using a Jeans 
Axisymmetric Multi-Gaussian mass model, we investigate the limit that a relaxed mass distribution 
produces all of the measured velocity dispersion, and find an upper limit on the black hole mass at 
2.0 ± 0.2 X 1O^M0. When assuming the young stars whose spectra we observe remain in a thin disk, 
we compare Keplerian velocity fields to the measured two-dimensional velocity field measured and 
fit for a mass profile containing a black hole point mass plus a radially- varying spherical component, 
which suggests a lower limit for the black hole mass of 8.7 ± 0.3 x lO^M©. Our measurements of the 
stellar velocity dispersion place this AGN within the scatter of the Mbh-<^* relation. As NGC 6240 is 
a merging system, this may indicate that the relation is preserved during a merger at least until the 
final coalescence of the two nuclei. 

Subject headings: Galaxies: kinematics and dynamics - galaxies: nuclei - galaxies: interactions - 
galaxies: individual (NGC 6240) 



1. INTRODUCTION 

Major mergers are thought to be an important factor 
in galaxy evolution. The scenario is as follows: when 
two gas-rich galaxies of comparable mass collide, large 
amounts of gas are tunneled into the central region, 
fuel ing active s tar formation and nu clear a ctivity (see 
e.g. iSanders fc M irabel 199 (1: iBarnes fc Her nauistjll996l: 
Genzel et allTToQSI: IDi Matteo et ahlEoOSt iHopkins et al.l 
20061) . During this phase of merging, the starburst- 
ing galaxy produces copious infrared emission from dust 
heated by young stars and by the active galactic nu- 
cleus (AGN). We see these galaxies in the local universe 
as (Ultra-)Luminous InfraRed Galaxies - (U)LIRGs. 
ULIRGs have infrared luminosities of more than 10^^ Lq. 
This burst of star-forming activity then uses up much of 
the gas; the remainder is blown out through a combina- 
tion of stellar winds and feedback from the AGN. The 
scenario then posits that after a major merger, gas-rich 
galaxies become gas-poor, star formation is largely ex- 
tinguished, and eventually a "red and dead" elliptical 
galaxy is produced with a more massive black hole at its 
core. 

NGC 6240 {z = 0.0243, d ^ 98 Mpc for Ha = 75 km 
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s-1 Mpc-i, 1" = 470 pc), with Ljr - IO^ ^Lq sits on 
the boundary between LIRGs and ULIRGs. Because of 
its close proximity and spectacular tidal tails and loops, 
it has become the prototypical example of a gas-rich sys- 
tem in the phase where the two nuclei are close to merg- 
ing into one. It has been studied in great detail and in 
almost every wavelength regime (e.g. x-ray - Komossa et 
al. 2003; optical - Gerssen et al. 2004; near-IR - Max et 
al. 2005, 2007; Scoville et al. 2000; Tecza et al. 2000; En- 
gel et al. 2010; mid-IR - Armus et al. 2006; mm - Tacconi 
et al. 1999; radio - Gallimore & Beswick 2004, Hagiwara 
et al. 2011). Near the core of NGC 6240, the nuclei of 
the two progenitors are visible 1.2-1.5 arcsec apart, de- 
pending on wavelength. Each of these nuclei holds an 
AGN; the two sources are res olved in hard x-rays by the 
Chandra X-Ray Observatory (jKomossa et al. 2003). The 
AGNs are deeply obscured at optical wavelengths, how- 
ever, due to large quantities of dust also present in this 
region . By looking into the near-infrared, iPoUack et al.l 
(j2007[) have seen young star clusters through some of the 
dust, products of the most recent close passage of the 
nuclei (also visible in Figure [T]). 

Supermassive black hole masses are known to scale 
with certain host gala xy properties, such as bulge 
light and mass (e .g. iK ormendv fc Richstond 119951 : 
iKormendv fc Gebha rdt 200l!: lMagorrian et al.lllQQSD and 
bulge stellar velocity dis persion ([Gebhardt et al. 020001 : 
iFerrarese fc Merrittll2000 i. Because these quantities can 
evolve significantly throughout the major merger pro- 
cess, black hole parameters may also evolve. It is not 
known whether the black hole mass grows as a result of 
the evolution in the host's bulge, or if the bulge growth 
is moderated by the process of AGN feedback; indeed, 
it is likely a combination of feeding and feedback pro- 
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Fig. 1. — Keck adaptive optics image of NGC 6240 in K' band 
(data first publisfied in lMax et al.li2007i ). Tfie dark-blue-enclosed 
regions in the north and south nuclei, which are separated by about 
1.6 arcsec, have each been re-scaled to highlight their interior struc- 
ture. The more diffuse image of the rest of the galaxy's nuclear 
region uses a logarithmic color map. Many i ndividual young sta r 
clusters can be seen exterior to the two nuclei llPoUack et al.ll2007l1 . 
In this image, north is up and east is to the left. 



cesses that maintains these scahng relations. In order to 
understand this coevolution, one needs to study systems 
that are currently merging and eventually to compare 
these observations to the most detailed merger simula- 
tions available. NGC 6240 represents the ideal candidate 
for such detailed observations. 

One important aspect of this effort for NGC 6240 is 
obtaining an accurate measurement of the black hole 
masses. A variety of methods are used to measure black 
hole mass, some of which we review here. One estab- 
lished method uses full three-integral modeling of stel- 
lar orbits; this metho d, most commonly used in ellipti- 
cal ga laxies, (see e.g. iSiopis et al.ll2009t iGiiltekin et al.l 
l2009al and references therein) uses line-of-sight velocity 
distributions and light profiles to create a detailed dy- 
namical mass profile of the galaxy, including the black 
hole mass. Two- and three-integral mo dels are also 
some t imes used on bulg e s of spirals (e.g. IDavies et al.l 
[2006t lOnken et al.l [2007t ICappellari et al.l 120091) . Re^ 
verberation mapping (e.g. iDennev et aL n2oo9r and ref- 
erences therein) is a method of measuring black hole 
masses in AGN; it uses time-resolved brightness fiuc- 
tuations in the continuum versus lines from the broad- 
line region to estimate the size of the broad-line re- 
gion and therefore the mass of the black hole power- 
ing the AGN. Masses measured with reverberation map- 
ping can be calibra ted to match the Mbh-i^* relation 
(lOnken et ar'200i). Bla ck hole masses have been de- 
rived ((Sh ields et a l] [20031 and references therein) from 
the continuum luminosity and the width of the broad 
H/3 line. The masses of supermassive black holes pow- 
ering quasars have been e stimated using x-ra y luminos- 
ity as an indicator, as in iKiuchi et all (|2009l ). In some 
cases, it is possible to measure the Keplerian rotation 



of ionized or molec ular gas around a b l ack hole and de- 
duce i ts mas s, as inlHarms et al.l (I1994D , iHicks fc MalkanI 
(|2008D . and iNeumaver et al.l ()2007D . In a few galaxies, 
masers in the disk allow precise velocity (and therefore 
black hole mass) measurements (Hcrrnst ein et al.) 119991 : 
IMivoshi"eral.. 1995). In the Galactic Center, individual 
stars have been resolved and their orbits around the black 
hole tracked ast rometrically over time to determine th e 
black hole mass (IGhez et al.ll2008l: IGillessen et al.l[2009l) . 
Recently, adaptive optics have been used to resolve stel- 
lar dynamics inside the sphere of influence of black holes 
in ne a rby galaxies (e.g. Davies et al. '2006'; 'Nowak _et al.l 
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technique which we now build upon. 

In a system such as NGC 6240, we are limited in our 
choice of method; the black holes are obscured by dust 
and the general system dynamics are unrelaxed because 
of the ongoing r nerger. Keck laser guide star adaptive 
optic s (LGS AO. lWizinowich et al.|[2006Hvan Dam et al.l 
|2006[) enables us to address both of these challenges. By 
looking in the near-IR, we look through much of the dust 
obscuring the relevant kinematics. The high spatial res- 
olution afforded by the adaptive optics system allows us 
to focus on stellar dynamics within the sphere of influ- 
ence of the black hole, unconfounded by the unrelaxed 
dynamics of the system at large. 

2. OBSERVATIONS AND DATA 

We began by observing NGC 6240 with the W. M. 
Keck II 10-meter telescope using the Near InfraRed Cam- 
era 2 (NIRC2, Pis - K. Matthews & T. Soifer) and the 
Keck LGS AO system to obtain high-resolution imaging. 
Our images, taken in the K' filter and using the narrow 
camera (wit h a 0.01" /p i xel pl ate s cale), were previously 
published in lMax et all ()2007[ ) and lPollack~al . (2007). 

We then observed NGC 6240 with the O H-Suppressing 
Infra Red Imaging Spectrograph (OSIRIS, iLarkin et al.l 
I2006D . on the W. M. Keck II telescope using LGS AO. 
OSIRIS is a near-infrared integral field spectrograph with 
a lenslet array capable of producing up to 3000 spec- 
tra at once. The spectral resolution ranges from about 
3400 in the largest pixel scale to 3800 in the three finer 
pixel scales; this resolution is sufficient to resolve spec- 
tral regions between the OH emission lines from Earth's 
atmosphere. Our data are comprised of two 600-second 
exposures in the Kn5 filter (2.292 /im - 2.408 /im) with 
the 0.035" /pixel plate scale taken on 21 April 2007. With 
this filter, we observe the CO (2-0) and (3-1) bandheads 
at 2.293 /^m and 2.323 /xm rest wavelength (2.345 /xm 
and 2.380 /im observed) respectively. Figure [5] shows an 
example spectrum. 

The Keck LGS AO system uses a pulsed laser tuned 
to the 589 nm Sodium D2 transition, exciting atoms in 
the sodium layer of the atmosphere (at ~95 km) and 
causing spontaneous emission. Thus, the laser creates a 
spot in the upper atmosphere which allows the AO sys- 
tem to monitor turbulence below the sodium layer via a 
Shack-Hartmann wavefront sensor and correct for it with 
a deformable mirror. A laser guide star enables high- 
order corrections to the wavefront, but relies on a nat- 
ural guide star (which may be fainter and farther away 
from the target than if no laser were used) to make cor- 
rections to image motion (tip and tilt). Our tip-tilt star 
(R=13.5 mag) is 35 arcseconds to the northeast of the 
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Fig. 2. — An example of the CO absorption bandhcads observed 
in the vicinity of the southern black hole with OSIRIS. This spec- 
trum was created by binning the light from a 3x3 pixel region 
around the black hole. The thick red line overplotted shows the 
results from template fitting to obtain kinematics, as detailed in 
Section [XT] The red dots along the horizontal line show the resid- 
uals of the fit about zero, shifted upwards to fit on the same plot. 

nuclei. To estimate the point-spread function (PSF), we 
took short exposures of our tip-tilt star before and after 
our observations; these will be described in Section 3.2. 

Our OSIRIS data were reduced with the 
OSIRIS Data Reduction Pipeline v2.2 (available at 
[http://irlab.astro.ucla.edu/osiris), which includes mod- 
ules to subtract sky frames, adjust channel levels, 
remove crosstalk, identify glitches, clean cosmic rays, 
extract a spectrum for each spatial pixel, assemble 
the spectra into a data cube, correct for atmospheric 
dispersion, perform telluric corrections, and mosaic 
frames together. 

3. METHODS 

Using the Kn5 filter, we observed the CO (2-0) and (3- 
1) bandheads at 2.293 fj,m. and 2.323 /itm rest wavelength 
respectively. These molecular features come from the 
atmospheres of later-type giants and supergiants. Stel- 
lar kinematics are less likely to be disrupted by non- 
gravitational forces than gas motions, giving us a poten- 
tially more robust measurement of the black hole mass. 

3.1. Measuring the Kinematics 

We begin by creating a signal-to- noise (S/N) map 
using two methods. We first calculate the S/N theo- 
retically, by adding noise components in quadrature: 
photon noise from the source, photon noise from the 
sky, read noise and dark current from characterizations 
of the detector. This produces an upper limit to the 
S/N because other noise sources may be present as well. 



S/N = 



F-ft 



^{Ft + Fstsii-^fh + [RN^ + (i)2 + Dt)7ij,,, 

In this equation, F is the flux from the galaxy, Fs is 
the flux from the sky, 7 is the detector gain, t is the 
exposure time of galaxy frames, ts is the exposure time 
of sky frames, RN is the readnoise of the detector, D 
is the dark current, and npix is the number of spectral 



pixels. For a mor e thoroug h discussion of this equation, 
see Section 9.9 in iMcLeani (|200a . 

It is most common to report S/N as an average signal- 
to-noise ratio per spectral pixel or per resolution element. 
We adjust the above equation appropriately by dividing 
the numerator and each variance {not the square- rooted 
noise components) by rispecpix, the number of spectral 
pixels included in the region of interest. In our brightest 
spatial pixels, we find a S/N per pixel of ~ 40. 

To confirm our theoretical calculations, we calculate 
S/N empirically from the spectra. Bluewards of the CO 
(2-0) bandhead, we have a spectral region of the galaxy 
uncontaminated by lines. Fitting for this continuum 
level, we find the signal present in each spectrum. After 
subtracting off the continuum fit, the root mean square 
of the residuals gives us a representation of the noise per 
pixel. We find good agreement between our empirical 
and our theoretical estimates of signal-to-noise ratios. 

Next we bin our data wi th optimal Voronoi tesselati ons 
using code developed by iCappellari &: Copinl (|2003[ ) to 
improve the signal-to-noise ratio in our fainter regions. 
This algorithm calculates a set of bin centroids according 
to specific criteria on the topology, morphology, and uni- 
formity of S/N of the final bins. That is, starting from 
the unbinned pixel with the highest S/N, the algorithm 
will include a pixel in the bin if it is adjacent to the start- 
ing pixel, does not significantly reduce the "roundness" 
of the bin, and, if the bin's S/N is too low, brings the 
S/N of the final bin closer to the chosen S/N threshold. 
In this way, pixels which already have a S /N at or above 
the chosen threshold are not binned and therefore spatial 
resolution is not sacrificed unnecessarily. In lower S/N 
regions, pixels are binned just enough to provide mean- 
ingful measurements. The morphological requirements 
create bins that are most likely to share similar v eloc- 
ities and dispersions. Following lEngel et all (|20TO ). we 
choose a S/N (per spectral pixel) threshold is 20 for each 
bin of pixels; the optimal bins are shown in Figure [31 

Once we have binned our spectra appropriately, 
we use the Penalized Pi xel Fitting code from 
ICappellari fc EmsellemI (2004) to fit radial velocities to 
each spectrum. This method implements a maximum 
penalized likelihood approach for extracting stellar kine- 
matics from absorption-line spectra. The algorithm 
parametrically expands the line-of-sight velocity distri- 
bution as a Gauss-Hermite series and allows for the 
choice of a penalty against higher-order moments. This 
penalty will bias the fit against higher-order moments, 
so that the fit must be improved by a certain specified 
amount to include them. In this way, higher-order fits are 
possible where spectra have a high enough signal-to-noise 
ratio, but the fits will tend to simple Gaussians in the low 
signal-to- noise limit. In order to trust higher- order mo- 
ment s, a signal-to-noise ratio of > 50 (jCappellari et al.l 
120091 ) is usually required. Since binning our spectra up 
to a S/N of 50 would decrease our spatial resolution sign- 
ficantly, we choose to fit simple Gaussians. 

One key feature of this code is the option of including a 
set of stellar templates from which to fit the kinematics. 
There has been some debate on which templates most ac- 
curately rep r esent the stellar populations of NGC 6240. 
iTecza et all ()2000[ ) conclude that the bulk of the light 
in the near-infrare d is due to lat e K or early M su- 
pergiants, while .Engel et al.l (|2010l ) argue that late-type 
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Fig. 3. — Top panel: The Voronoi tesselation bins we impose on 
our data to equalize to S/N of 20 in the spectrum associated with 
each spatial bin. The bins near the center are small (one pixel) 
because the flux is high enough that the S/N is already above the 
threshold. In the outer regions, pixels have been binned together 
so that the spectra have a S/N of ^20. Bottom panel; Another 
representation of the S/N of bins. Each point represents a bin; 
the plus symbols are bins that contain only one of the original 
pixels — they lie above the threshold line because they had sufficient 
S/N initially. The open squares represent bins whose members 
originally had insufficient S/N; they were binned together until 
they approximately reach the threshold value. 

giants characterize the global stellar population better. 
With Cappellari's method, we are able to input a vari- 
ety of stellar templates and allow the parametric fit to 
select the best combination of templates for each spec- 
trum. In our final iteration of the code, we selected 5 
stella r templates from the GNIRS library (IWinge et al.l 
l2009f) that exhibited the deepest CO bandheads. The 
stellar templates chosen were: HDI12300 (spectral type 
M3III), HD198700 (Kill), HD63425B (K7III), HD720, 
(K5III), and HD9138 (K4III). The region of these spec- 
tra around the CO bandheads is shown in Fig HI 

We show the measured velocity and velocity dispersion 
maps from OSIRIS in Figures [S] and IH respectively. We 
estimate the errors in the velocity and velocity disper- 
sion measurements by running a Monte Carlo simulation, 
adding the appropriate amount of random noise to each 
spectrum and refitting them 100 times. Formal mea- 
surement errors in our velocity and velocity dispersion 
measurements are ~ 10 km s~ in the brightest regions, 
and average about 30 and 70 km s^^, respectively, over 
the whole fitting region. 

With velocity and velocity dispersion maps, it is tempt- 
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Fig. 4. — Spectra of the five late- type giants and supergiants 
used as stellar velocity templates in our dynamical modeling. The 
deep CO bandheads match those seen in our spectra (Figure [2]|. 
The templates were shifted in wavelength to correct for peculiar 
velocities of the stars, to match the CO transitions in a vacuum. 

ing to think about the measured vja of the system, to 
determine how much of the stellar dynamical energy is 
in rotation. However, this must be handled carefully in 
a system like NGC 6240. In our data, we measure ve- 
locity peaks of ~ 200 km s~^, and velocity dispersion 
that varies between roughly 200 and 300 km s~^, sug- 
gesting that the energy is approximately evenly divided 
between ordered rotation and random orbits, with per- 
haps up to 50% more in the latter. However, in this 
region of NGC 6240, we know there are clumps of inter- 
veni ng material such a s spiral arm remnants or tidal tails 
(e.g. lEngel eraII[20Tol ) . which would cause an increase m 
the measured velocity dispersion that is not yet inherent 
to either the nuclear disk or the spheroid. Because we do 
not have a way of determining how much of the velocity 
dispersion measured is from relaxed material, our analy- 
sis yields only a lower bound on vja. We conclude that 
a significant fraction of the kinetic energy is in rotation. 

Once we have maps of velocity and velocity disper- 
sion, we compare our data to models that contain a black 
hole. Each model, of course, comes with its own set of 
assumptions, and in a late-stage merging system such 
as NGC 6240, we must think carefully about what such 
models can tell us. While the quality of these data is 
clear, understanding how to analyze them is not straight- 
forward. Here we discuss two possible ways to measure 
the black hole mass from these data and compare the 
results of the two methods. 

3.2. Dynamical Analysis - JAM Modeling 

We begin our analysis by utilizing the JAM modeling 
code (Jeans An isotropic Multi-G aussian expansion dy- 
namical models [CaDpellarill2008D . a technique based on 
the two- integral axisymmetric Jeans formalism but which 
has been expanded to allow for anisotropy via the param- 
eter /J^ = 1 — (cTz/fTfl)^. JAM modeling efficiently utilizes 
the axisymmetric dynamics seen near the south nucleus, 
and does not require higher-order Hermite moments. 

This method, which fits w,ms = vi? + o^, is likely to 
overestimate the black hole mass by assuming the dy- 
namics measured belong to a relaxed system. In fact. 
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Fig. 5. — Stellar velocity field measured from the OSIRIS IFU 
data near the southern black hole (red x). Typical errors arc ~ 10 
km s~^ in the brightest regions and ~ 30 km s~^ overall. Large 
pixels on the periphery have been binned using Voronoi tesselation 
to improve S/N (see Figure[3ll. Voronoi bins with centroids further 
than ~ 115 pc from the black hole (in projection) have been colored 
white to mask out regions which are less affected by the black 
hole's gravity. (The data used to create this figure are available for 
download in the online journal.) 
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Fig. 6. — Stellar velocity dispersion measured from the OSIRIS 
IFU data near the southern black hole (red x) . Typical errors are ~ 
10 km s~^ in the brightest regions and ~ 70 km overall. Large 
pixels on the periphery have been binned using Voronoi tesselation 
to improve S/N (see Figure [Sjl. Regions far from the black hole 
have been masked out, as in Figure [5] (The data used to create 
this figure are available for download in the online journal.) 

an unknown fraction of the measured velocity disper- 
sion is due to intervening material, such as tidal tails, 
that has not yet reached dynamical equilibrium. Addi- 
tionally, this method assumes axisymmetric and smooth 
light and mass profiles. Finally, we assume a constant 
mass-to-light ratio. 

The JAM modeling code requires a high-resolution 
light profile, which we parametri ze using the M ulti- 
Gaussian Expansion code (MGE Cappellaril l2002f ) de- 
signed to work with the JAM code. To fit our light 



profile over a larger field of view than is available in our 
OSIRIS data, we use our K' NIRC2 imaging and mirror 
it about the minor axis of the nucleus. We do this be- 
cause the southeast side of the nucleus is considerably 
less extincted, so our signal-to-noise ratio is much im- 
proved. Once we have symmetrized the observed light 
profile, we the n de-extinct this u sing the extinction map 
of Figure 9 in lEngel et al.l ()201Cl[ ). Using this extinction 
map, our recovered intrinsic brightness peaks at the same 
location as the kinematic center of our dynamical data, 
which means that in K-band the nucleus is only partially 
extincted. (In contrast, the extinction at visible wave- 
lengths is so severe that the entire regi on surrounding t he 
south black hole cannot even be seen (|Max et al.ll2"005l) .') 
This is important, as most of the orbital information re- 
ported by JAM modeling is contained in the light profile. 
In order to determine the appropriate mass-to- light ratio, 
anisotropy parameter /3z and black hole mass, we com- 
pare the resulting dynamical models to those measured 
from our OSIRIS data in K-band. 

Our best-fit models for this method measure a black 
hole mass of 2.0±0.2 x IO^Mq, and are shown, along with 
the symmetrized Vrms data for comparison, in Figure [71 
The reduced statistic, fitting over 120 points, is 2.21. 
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Fig. 7. — Left: Symmetrized Vrms map (as in 'Cappellaril'20Q3) 
from OSIRIS observations. Center: Best-fitting axisymmetric 
JAM model, on the same scale and color bar as data (left). Right: 
Map of residuals (absolute value of difference between left and cen- 
ter panels) with new color bar to the (right). In all panels, the black 
contours show the Gaussian expansion of the light profile, and the 
black dots represent the centroid of each spatially-binned spectrum, 
where velocity and velocity dispersion are measured. Axes show 
distance from the central black hole in arcseconds. 



3.3. Dynamical Analysis - Thin Disk 

As a sanity check on the JAM model, and to provide 
a lower-limit to the black hole mass, we explore a simple 
model comparing the velocity field to that of a thin disk 
with a given enclosed mass profile exhibiting Keplerian 
rotation. 

In this method we do not include a dispersion compo- 
nent because intervening material may inflate that mea- 
surement. Here we assume that the energy in intrinsic 
dispersion in the nuclear stellar disk is negligible com- 
pared to the energy in rotation. This should give a lower 
limit to the black hole mass. 

We are measuring the dynamics of young stars in 
the very nucleus of a gas-rich merger. A thin disk 
of young stars may be expected because they can 
form out of the nuclear disks of g as and dust seen i n 
merger si mulation s by, fo r example, Mayer et al." ('20071); 
iKazantzidis et al.l (|2005[) : iHopkins fc Ouataert (,2010al) 
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and i n observations by iRiffel fc Storchi-Bergmamil 
(|20T0l) : iHoHkins fc Quatae^ (|2010b[) . We clearly see in 
the OSIRIS velocity field a sharp steepening of the ve- 
locity gradient in the region of the black hole: a sign of 
strong rotation, and reminiscent of a thin Keplerian disk 
embedded in a larger spherical potential. 

We focus our mass modeling on the region of the sphere 
of influence of the black hole, where the dynamics are 
most likely to be well-behaved. We begin with a thin 
Keplerian disk model, using v = ^JGMbh /t- This will 
attribute all mass enclosed to the black hole, and there- 
fore overestimate the mass; we only use this case to check 
the sub-pixel position of the black hole and get a rough 
set of parameters over which to fit a more complex model, 
which includes a spheroid, described in Section [3.3.11 

Once a model velocity field is constructed for a spe- 
cific set of parameters, we use that information to create 
a synthetic datacube: for each spaxel, we use a template 
CO bandhead spectrum shifted in velocity space to the 
appropriate velocity and weighted by the total flux in 
that spaxel. With this synthetic datacube we perform 
a wavelength-by-wavelength convolution with the PSF 
to simulate the residual smearing from the atmosphere 
and optical system. Convolving each wavelength slice 
with the PSF models the true image blurring, as op- 
posed to smoothing the overall velocity field. Once the 
datacube has been convolved in this way, we remeasure 
the velocities from each spatial pixel's spectrum to derive 
the smoothed velocity field. We compare the resulting 
velocity field with the observed velocities measured by 
OSIRIS. We repeat this process, making models for vari- 
ous parameter sets, and compare each with the observed 
data. 

3.3.1. Model Parameters 

At these spatial scales, the velocity field of a thin Ke- 
plerian disk around the black hole can be fully described 
by a few parameters: black hole position, disk inclina- 
tion, disk position angle, and enclosed mass. We fit 
the position angle of the velocity field first, as it may 
be fit largely independently of the other black hole pa- 
rameters. Our tests indicate that the positive velocity 
peak falls at a PA of 130° measured counter-clockwise 
from north, relative to the black hole's position. We 
adopt this position angle in our later fits. The position 
of the black hole is not trivial to pinpoint. Because of 
6-8 magnitu des of extinction di stributed unevenly across 
the nucleus (|Engel et al.ll2010D . we avoid using our K'- 
band isophotes to find the center of the southern nucleus, 
and instead use the measured velocities to determine the 
kinematic center. We have a close estimate of the posi- 
tion of the southern black hole using relative astrometry 
from the northern black hole, whose position is visible in 
our data ((Max et al..,2007i) . To improve our positional 
accuracy, we allow the black hole position to vary in the 
fit by sub-pixel amounts in a simple code that only takes 
into account the black hole mass. We then fix the po- 
sition of the black hole, along with our PA, for our full 
suite of models. 

As AO measurements can be affected by the PSF 
achieved in our observations, we must carefully consider 
the PSF used in our models. We use short exposures 
of our tip-tilt star bracketing our observations which are 
useful to characterize the performance of the AO system 



and the conditions during the evening. Our first PSF 
model is a Moffat fit of the tip-tilt star, which has a Strehl 
ratio of 20% and FWHM of 65 mas. We find that the 
north nucleus point source has a FWHM of 63 mas, con- 
sistent with that of the tip-tilt star. However, the tip-tilt 
star gives a lower-limit to the actual Strehl ratio since, 
for shorter exposure times appropriate to the tip-tilt star, 
the l ow-bandwidth wavefront sensor (jWizinowich et aLl 
|2006[ ) did not have time to settle into the most accurate 
correction. The longer exposures on NGC 6240 itself do 
provide sufficient settling time. This effect is partly offset 
by additional blurring due to the offset between the tip- 
tilt star and NGC 6240 (anisokineticism) . Because the 
distance between the tip-tilt star and th e nuclei is only 
35" , this effect is expected to be modest ()van Dam et al.l 
12006( 1. and does not completely counteract the aforemen- 
tioned effect from the low-bandwidth wavefront sensor. 
We see that the AO system performs better during our 
longer galaxy exposures. To accommodate this, we ad- 
just the Moffat fit coeffecients to Strehls of 25%, 30%, 
35% and 40%, while maintaining the FWHM (see radial 
profiles in Figure [8]). The Strehl ratio then becomes a 
parameter in our model-fitting procedure, in which we 
create models with each of these PSFs to compare with 
our data. 



PSF 40% Sir 
PSF 35% Str 
PSF 30% Str 
PSF 25% Str 
PSF 20% Str 




0.4 

Radius in arcsec 

Fig. 8. — The five PSFs used by our model fitting routine. The 
Strehl ratio is varied as a free parameter in order to avoid sys- 
tematic errors in black hole mass measurement based on PSF 
mismatch, while maintaining a FWHM of 65 mas, to match the 
FWHM measured in the data. 



In an asymmetric, dusty system such as NGC 6240, 
we must also carefully consider the region over which we 
compare our models to the data, masking out regions 
that are less important to match. We expect that our 
model will most accurately describe the velocities close to 
the black hole. At some distance from the black hole, we 
expect other galaxy and tidal components to dominate, 
which may affect the model fit. To test this concern, we 
use a variety of different masks to vary the distance from 
the black hole at which we stop our comparison. We find 
that the best-fitting enclosed mass is not sensitive to the 
mask selection. 

Our mass profile includes a radially-varying contribu- 
tion to the mass profile, representing a spheroidal com- 
ponent of the galaxy. We construct models that in- 
clude both a point mass at the center and a spherically- 
symmetric mass profile to mimic the inner regions of 
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the bulge. We parametrize the spheroidal component 
as p{r) = por^ , where p is the mass density. We fit three 
parameters to the enclosed mass profile: the black hole 
mass Mbh, the normalization for the spheroidal com- 
pontent pq, and the power-law index of the spheroidal 
component, 7. 




-50 50 100 -50 50 100 

Distance from BH (pc) Distance from BH {pc} 



Fig. 9. — (a) The observed velocity field of our OSIRIS data for 
comparison with our best-fit models, (b) The best-fitting velocity 
field for a mass model containing a point mass and an extended, 
spherically-symmetric component, (c) The residuals found by sub- 
tracting the observed velocity field from the model showrn in (b). 
Note that the model yields a good fit to the central steep velocity 
gradient associated with the black hole. In each panel regions far 
from the black hole have been masked out, as in Figure [5] 



3.3.2. Results of Thm Disk Model Fitting 

Our best fit gives a black hole mass of 8.7x IO^Mq with 
a spheroid of 2.9 x 10® Mq within 100 pc of the black hole 
and an index of 7 = 1.5 (see velocity field of this model 
in Figure El). The reduced statistic for this model is 
3.7, fitting over 65 velocity datapoints. To estimate the 
accuracy of this measurement, we perform a Monte Carlo 
simulation, fitting models to our observations with 100 
different representations of noise added. The distribution 
of fitted masses can be well-fit by a Gaussian, the width 
of which represents our one-sigma error bars. Including 
this, our best fit shows a system with a black hole mass 
8.7 ±0.3 X 10^ Mq. 

3.3.3. Trends in the Reduced Map 

It is also instructive to look at trends in the reduced 
maps. The simplest example is shown in Figure IIOI 
demonstrating that, for a thin disk model such as this 
one, the black hole must be more massive to match the 
velocity field when the inclination is lower, for i < 60°. 
This map also shows that the black hole mass increases 
again at the highest inclinations, as the model tries to 
match the velocity field not only along the major axis but 
in the surrounding spatial pixels as well (the well-known 
"spider diagram" ). If NGC 6240 were dust- free, the mor- 
phology of the region around the south nucleus could, 
in principle, be used to independently measure (or con- 
strain) the inclination. However, the northwest half of 



the circumnuclear disk is heavily obscured in NGC 6240. 
In this case, the image can only provide a lower-limit for 
the inclination 55°). 




Inclination (°) 

Fig. 10. — The map of reduced statistics for our thin disk 
model as the black hole mass (vertical axis) and the inclination 
(horizontal axis) are varied. The remaining parameters (spheroid 
paramaters po and 7, velocity offset, and PSF) are held fixed. We 
see that as the inclination decreases (the disk of stars becomes more 
face-on) , a larger black hole mass is required to match the observed 
radial velocity peaks for inclinations < 60°. As the inclination 
increases beyond 60° (the disk of stars becomes more edge-on) , we 
also see the best-fit black hole mass increase; this is likely because 
a larger black hole is required to affect off-axis regions of the spider 
diagram when the disk is more edge-on. These fits were measured 
with 65 datapoints. 

We see a demonstration of the trade-off of mass be- 
tween the spheroid and black hole by looking at a 2- 
dimensional map of reduced statistics (Figure llip . 
varying the amount of mass that goes into each compo- 
nent of our mass profile. As more mass is put in the black 
hole, models with less mass in the radial component fit 
better; as less mass is put in the black hole, the mass in 
the radial component must increase to best fit the data. 
In each case, the combined mass reaches a maximum at 
~ 2 X 1O^M0 within 100 pc of the black hole. 

Lastly, we consider our choice of varying PSFs in our 
model fitting. Our best-fit models from this technique 
prefer PSFs with a Strehl ratio of 25%, near the lower 
limit of what we would expect from our data. We com- 
pare these to models with 40% Strehl, approximately the 
best PSF we could hope for under the conditions of these 
observations. When forcing an improved PSF, the best- 
fitting black hole mass increases slightly, to 9.2 x IO^Mq. 
This is only slightly outside of the 1-a range, and has a 
reduced statistic of 4.8, considerably higher than our 
better-fitting models. While we are confident that fitting 
for PSF was a good choice, we also note that the black 
hole mass is not very sensitive to this parameter. 

4. DISCUSSION 

4.1. How Much Mass Could Be Due to a Nuclear Star 

Cluster? 

It is important to consider the limitations of our ap- 
proaches. With a black-hole-only thin disk model, one 
does not directly measure the black hole mass; one mea- 
sures the mass enclosed in the central OSIRIS pixel, 17 
parsecs on a side. Including a spheroidal component im- 
plicitly assumes that the stellar mass density profile is 
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Fig. 11. — The map of reduced statistics as tlie black hole 
mass component (horizontal axis) and the spherically-symmetric 
mass component (vertical axis) are varied. The other parameters 
(inclination, density profile index 7, and velocity offset) are held 
fixed. We see that on the right side (at high black hole mass), the 
best fitting models have less mass in the spheroidal component. On 
the left side (with black hole mass becoming negligible), the mass 
of the spheroidal component flattens out at 2 X 10® Mq, the total 
enclosed mass at our resolution limit. These fits were measured 
with 65 datapoints. 

smooth and that ah other mass is due to the black hole. 
What fraction of this mass might be due instead to a 
nuclear star cluster? 

To address this concern, we refer to our high-resolution 
NIRC2 K' imaging. Through a careful deconvolution by 
the PSF, we arrived at a cleaned image of the south nu- 
cleus. We flux-calibrated this image by match ing the 
large-scale luminosity to that reported in lEngel et al.l 
(j201C|. and assumed their mass-to- light ratio of 1.9. 

This analysis showed that there could be up to 3 x 
lO^Msun of stellar mass within the sphere of influence 
of the black hole, consistent with the mass in our fitted 
spheroidal component. 

4.2. The Mbh — cr* Relation 

With our measurement of the black hole mass, it is 
interesting to consider where NGC 6240 would fall on 
the Mbh — relation, which compares the black hole 
mass to the stellar velocity dispersion of the bulge of the 
host galaxy. Because NGC 6240 is a merging system, it 
plausibly lacks a relaxed bulge component; therefore it 
is not easy to define exactly where or how the velocity 
dispersion should be measured to compare to the Mbh — 
tr* relation. Still, it is worth looking at merging active 
systems since they are in the process of evolving along 
these relations. Here we are able to take our first glimpse 
of where a system might fall on these relations while in 
the process of merging. Does the black hole grow more 
quickly than the larger-scale galactic properties, or must 
it play catch-up after the galaxy's bulge has settled back 
to an equilibrium state? 

The Mbh — cr* relation is well-defined only in systems 
that are dynamically-relaxed on the large scale. In such 
a system, the central black hole mass can be compared 
to the integrated velocity dispersion of stars in the bulge 
within one effective radius. NGC 6240 has two black 
holes; here we only consider the southern one. It is not 
trivial to measure the equivalent stellar velocity disper- 
sion associated only with this black hole in the southern 
galaxy, as the bulges of both progenitor galaxies have be- 



gun to merge. We expect the stellar velocity dispersion 
at this stage to be low compared to the final value; quite 
a bit of energy is still in ordered rot ation an d has y et to 
be randomized. As reported in Eng el et al.l ()2010t ). the 
south nucleus shows a maximum rotational velocity of 
^ 300 km s^^, putting as much or more energy in or- 
dered rotation as in dispersion. Our spectra also give us 
only the dynamics of l a.ter-type giants and s uperg iants, 
via the CO bandheads. iRothberg fc Fischei) (|2010l ) have 
suggested that this also gives an underestimate of veloc- 
ity dispersion as compared to measurements made using 
the Ca triplet absorption lines at 8500 A. 

Our spectra also only give us dynamics very close to 
the nucleus. We bin our OSIRIS data into one spectrum 
encompassing our entire south nucleus (inside a radius 
of ~ 300 pc); the measured velocity dispersion from the 
CO bandheads is 282 ± 20 km s~^. For comparison, we 
make the same m easurement using SINFONI data from 
lEngel et"al] (|20100 : this velocity dispersion, 310 ± 12 km 
s"-'^, encompasses a wider region (r < 500 pc). We also 
compare these two values of stellar velocity dispersion to 
previo us measurements from the literature. iTecza et al.l 
(|2000D measured ti* ~ 236 km s~^ using th e CO band- 
heads within 235 pc of the south nucleus. lOliva et al.l 
()1999t ) measured the stellar velocity dispersion integrated 
over the entire system, but used three different absorp- 
tion lines: Si 1.59 fim (313 km s^i), CO 1.62 fim (298 
km s~^), and CO 2.29 /im (288 km s~^). Estimates of cr* 
made from integrating over a larger fraction of the galaxy 
could be higher because they combine multiple dynam- 
ical populations (e.g. the north and south nuclei and 
intervening spiral arms); however, such estimates do not 
distinguish between material associated with this black 
hole and the northern black hole, as both are partially 
within the seeing-limited PSF. 

We plot our black hole mass measurement on the 
Mbh — o'* relation along with other dynamical black hole 
mass measurements compiled by Giiltekin et al. (2009l|) 
in Figure [12] We plot separate points for our black 
hole mass measurements under two different assump- 
tions. Because there is some ambiguity about the ap- 
propriate way to measure the stellar velocity dispersion 
in the bulge in a system such as NGC 6240, we plot our 
measured black hole mass with several different ct* values 
from the literature as well as with our measured cr* very 
close to the black hole. 

NGC 6240 appears to lie well within the scatter of the 
Mbh — cr* relation, which may suggest that the black 
hole mass and the bulge velocity dispersion grow simul- 
taneously and at similar rates during a major merger. 
NGC 6240 is a late-stage merger, but these data could 
also suggest that a system doesn't evolve along the 
Mbh — ^* relation until the very end stages of a merger, 
perhaps during nuclear coalescense. To test this, more 
systems at this and later stages of merging will have to 
be studied. This is interesting t o compare to galaxy 
merger simulations; for example, iDasvra et al.l ()2006f ) 
conclude that, if the accretion efficiency stays constant, 
the Mbh — relation should be maintained from mid- 
way between the first encounter and coalescence through 
to final relaxation. New higher-resolution simulations are 
being performed by Stickley et al. (in prep) which show 
that the stellar velocity dispersion stays low in the nu- 
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F ig. 12. — Plot of the M bh — relation as recently recalculated 
bv lGiiltekin et aLl Il2009bl ). The small black points and error bars 
represent the black hole mass measurements from the literature. 
The larger colored squares and diamonds represent our black hole 
mass measurements in the south nucleus of NGC 6240, using differ- 
ent values of cr* . Squares show the black hole mass measurement 
from JAM modeling, and represent an upper limit to the black 
hole mass. Diamonds show the black hole mass measurement from 
our thin disk approximation, and represent a lower limit to the 
black hole mass. The dark green points are our mass mea sure- 
ment paired with measurements of cr« from lTecza et all 1120001 ) and 
lOliva et al.l II1999I ). The blue and brown points plot our measure- 
ments of cr« from OSIRIS and SINFONI respectively. See text for 
further details. 



clear regions much longer than in the rest of the system; 
results will follow young and old stellar populations sep- 
arately to help understand the dynamics measured from 
specific lines. 

5. CONCLUSIONS 

Black holes are an important ingredient in galaxy evo- 
lution, as seen by the tight correlations between black 
hole mass and host galaxy properties. Hydrodynamic 
simulations of galaxy mergers have suggested that gas- 
rich mergers, which can drive galaxy evolution, can pro- 
vide fuel for black hole accretion. In turn, these accreting 
black holes can radiate enough energy to affect the sur- 
rounding galaxy. In light of these discussions, it is partic- 
ularly interesting to study systems that appear to be in 
the middle of such an evolutionary event. NGC 6240 rep- 
resents a nearby galaxy such as this: a merging gas-rich 
system with two actively accreting black holes. However, 
such systems are notoriously difficult to study because of 
their unrelaxed dynamics and their dusty cores. 

We have presented high-spatial resolution kinematics 
within the sphere of influence of the black hole in the 
south nucleus of NGC 6240, a nearby late-stage merger, 
made possible by laser guide star adaptive optics. For 
this test case, we have explored two possible methods for 
measuring black hole mass in such a system and com- 
pared their results and assumptions. 



We ha ve utilized the JA M modeling technique made 
public bv lCappellaril (|2008[ ). demonstrating that it is pos- 
sible to complete such an analysis on a late-stage merger. 
We point out that this technique likely overestimates the 
black hole mass by assuming that all measured velocity 
dispersion is due to a relaxed system (ignoring interven- 
ing unrelaxed material, e.g. tidal tails), and therefore 
report an upper limit of 2.0 ± 0.2 x IQt'Mq. To pro- 
vide a lower limit to the black hole mass, we explore 
the opposite assumption: that all velocity dispersion is 
caused by intervening material, and that the young stars 
sit in a thin disk around the black hole. This model 
of Keplerian rotation around a black hole plus smooth 
spheroidal mass profile suggests a black hole which is at 
least 8.7 ±0.3 x IO^Mq. 

We find that these two techniques provide measure- 
ments that are roughly consistent, and that follow the 
biases implied by their intrinsic assumptions. To deter- 
mine which set of assumptions is more reliable would 
require a detailed study of high-resolution galaxy simu- 
lations, which is beyond the scope of this paper. Still 
we are encouraged that, to within a factor of about two, 
both measurements agree. 

While we cannot make generalizations on how all black 
hole-galaxy coevolution must proceed, it seems that in 
this case, the black hole and the host galaxy parameters 
grow together along the Mbh — scaling relation, in- 
stead of one preceeding the other. We are beginning an 
observing campaign to study other local merging galaxies 
using the same techniques, to investigate a larger sample 
size. 
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